Competition of electron capture and beta-decay rates in supernova collapse 
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We calculate supernova electron capture and decay rates for various p/-shell nuclei using large- 
scale shell model techniques. We show that the centroid of the Gamow- Teller strength distribution 
has been systematically misplaced in previous rate estimates. Our total electron capture rates are 
significantly smaller than currently adopted in core collapse calculations, while the total f3 decay 
rates change less. Our calculation shows that for electron-to-baryon ratios Y e = 0.42-0.46 /3 decay 
rates are larger than the competing electron capture rates. 



PACS numbers: 26.50.+X, 23.40.-s, 21.60Cs, 21.60Ka 



I. INTRODUCTION 

Weak interaction processes play a decisive role in the 
early stage of the core collapse of a massive star (jl^). 
First, electron capture on nuclei in the iron mass re- 
gion, starting after the core mass exceeds the appropriate 
Chandrasekhar mass limit, reduces the electron pressure, 
thus accelerating the collapse, and lowers the electron-to- 
baryon ratio, Ye, thus shifting the distribution of nuclei 
present in the core to more neutron-rich material. Sec- 
ond, many of the nuclei present can also f3 decay. While 
this process is quite unimportant compared to electron 
capture for initial Y e values around 0.5, it becomes in- 
creasingly competative for neutron-rich nuclei due to an 
increase in phase space related to larger Qp values. How- 
ever, (3 decay on nuclei with masses A > 60 have not yet 
been considered in core collapse studies [[|. This is sur- 
prising since Gerry Brown pointed out nearly a decade 
ago M that certain nuclei heavier than A — 60, like 63 Co 
and Co, have very strong /3-decay matrix elements mak- 
ing it conceivable that they can actually compete with 
electron capture. Brown argued that this might have 
quite interesting consequences for the collapse. During 
this early stage of the collapse, neutrinos produced in 
both electron capture and (3 decay still leave the star. 
Therefore, a strong /3-decay rate will cool the star with- 
out lowering the Y e value. As a consequence, the Y e value 
at the formation of the homologous core (after neutrino 
trapping) might be larger than assumed. This results in 
a smaller envelope, and less energy is required for the 
shock to travel through the material. 

Following Brown's suggestion, /3 decay rates for nuclei 
in the mass range A — 48 — 70 were investigated ||. 
These studies were based on the same strategy and for- 
malism as already employed by the pioneering work in 
this field by Fuller, Fowler and Newman (commonly ab- 
breviated by FFN) The important idea in FFN was 
to recognize the role played by the Gamow- Teller reso- 
nance in (3 decay. Other than in the laboratory, (3 decay 
rates under stellar conditions are significantly increased 



due to thermal population of the Gamow- Teller back res- 
onance in the parent nucleus (the GT back resonance are 
the states reached by the strong GT transitions in the in- 
verse process (electron capture) built on the ground and 
excited states, see |^,||) allowing for a transition with a 
large nuclear matrix element and increased phase space. 
Indeed, Fuller et al. concluded that the (3 decay rates un- 
der collapse conditions is dominated by the decay of the 
back resonance. In a more recent work, Aufderheide et al. 
came to the same conclusion. Inspired by the indepen- 
dent particle model, the authors of Rcf. |3| estimated the 
f3 decay rates in a similar fashion as the e ectron capture 
rates, and phenomcnologically parametrized the position 
and the strength of the back resonance. This estimate 
was supplemented by an empirical contribution, placed 
at zero excitation energy, which simulates low-lying tran- 
sition strength missed by the GT resonance. 

When extending the FFN rates to nuclei with A > 60, 
Aufderheide et al. found indeed that the j3 decay rates 
are strong enough to balance the electron capture rates 
for Y e k 0.42 — 0.46. Nevertheless, these results have 
never been explored in details in core collapse calcula- 
tions. 

In recent years the parametrization of the electron cap- 
ture and (3 decay rates as adopted by FFN and Aufder- 
heide et al. have become questionable due to experi- 
mental data (7 -11 1 and have been critizised on the basis 
of more elaborate theoretical models ]l^-|l5|]. We will 
show in this paper, that although the previous weak in- 
teraction rates for core collapse are systematically incor- 
rect, the important observation that electron and (3 de- 
cay rates balance each other for a certain range of Y e 
values is indeed correct. Our conclusions will be based 
on large-scale shell model calculations for several key nu- 
clei which, due to Aufderheide et al. ||, contribute most 
significantly to the electron capture and (3 decay rates at 
various stages of the collapse. These shell model calcula- 
tions reproduce the measured GT strength distributions 
for nuclei in the mass range A = 50 — 64 very well [|l6| . 
Furthermore modern large-scale shell model calculations 
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also agree with measured half-lifes very well. Thus for 
the hrst time one has a tool in hand which allows for a 
reliable calculation of presupernova electron capture and 
(3 decay rates. Modern shell model calculations come in 
two varieties: large-scale diagonalization approaches fl7| 1 
and shell model Monte Carlo (SMMC) techniques 
The latter can treat the even larger model spaces, but has 
limitations in its applicability to odd-A and odd-odd nu- 
clei at low temperatures, which does not apply to the 
former. More importantly the diagonalization approach 
allows for detailed spectroscopy, while the SMMC model 
yields only an "averaged" GT strength distribution which 
introduces some inaccuracy into the calculation of the 
capture and decay rates. 

We will consistenly use in the following the shell model 
diagonalization approach to study these rates. Due to 
the very large m-scheme dimensions involved, the GT 
strength distributions have been calculated in truncated 
model spaces which fulfill the Ikeda sum rule. However, 
at the chosen level of truncation involving typically 10 
million configurations or more, the GT strength distri- 
bution is virtually converged. As residual interaction we 
adopted the recently modified version of the KB3 interac- 
tion which corrects the slight inefficiencies in the KB3 in- 
teraction around the N = 28 subshell closure [^0| . In fact 
the modified KB3 interaction i) reproduces all measured 
GT strength distributions very well and ii) describes the 
experimental level spectrum of the nuclei studied here 
quite accurately (l^^^l ■ As Otko shell model calculations, 
i.e. calculations performed in one major shell, overesti- 
mate the experimental GT strength by a universal factor 
pi] [23[ , we have scaled our GT strength distribution by 
this factor, (0.74) 2 . 

Large-scale shell model calculations of the electron cap- 
ture rates for key nuclei in the presupernova collapse 
have been reported already elsewhere |l4|]l5| ] ; see also the 
SMMC results in |b|] . In these studies it became appar- 
ent that the phenomenological parametrization of the GT 
contribution to the electron capture and (3 decay rates, as 
introduced in FFN || and subsequently used by Aufder- 
heide et al. |J, is systematically incorrect. These authors 
have placed the centroid of the GT strength distribution 
at too high an excitation energy in the daughter nucleus 
for electron capture on even-even nuclei, while for cap- 
ture on odd-A and odd-odd nuclei they underestimated 
the energy of the GT centroid noticeably. For capture on 
even-even nuclei this has comparably little effect as FFN 
overcompensate the misplacement of the GT centroid by 
a too large empirical contribution at zero excitation en- 
ergy. Typically the FFN rates are roughly a factor of 5 
larger than the shell model rates for capture on even-even 
nuclei like 56 > 58 Ni or 58 Fe. For capture on odd-odd and 
odd-A nuclei the misplacement of the GT centroid makes 
the FFN rates about 1-2 orders of magnitude too large 
compared to the shell model rates. As a consequence, 
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FFN and Aufderheide et al. have noticeably overesti- 
mated the electron capture rates for the early stage of 
the supernova collapse. 

Which consequences do the misplacement of the GT 
centroids have for the competing (3 decays? In odd-A and 
even-even nuclei (the daughters of electron capture on 
odd-odd nuclei), experimental data and shell model stud- 
ies place the back-resonance at higher excitation energies 
than assumed by FFN and Aufderheide et al. Corre- 
spondingly, its population becomes less likely at the tem- 
peratures available during the early stage of the collapse 
(Tg w 5, where Tg measures the temperature in 10 9 K) 
and hence the contribution of the back-resonance to the 
(3 decay rates for even-even and odd-A nuclei decreases. 
Due to Aufderheide et al. [^), some of the most impor- 
tant (3 decay nuclei (defined by the product of abun- 
dance and (3 decay rate and listed in Tables 18-22 in (|]) 
are odd-odd nuclei. For these nuclei, all available data, 
stemming from (n,p) reaction cross section measurements 
on even-even nuclei like 54 < 56 - 58 Fe or 58,60,62,64 Ni; and all 
shell model calculations indicate that the back-resonance 
resides actually at lower excitation energies than previ- 
ously parametrized. Consequently, the contribution of 
the back-resonance to the (3 decay rate of odd-odd par- 
ent nuclei should be larger than assumed in the compila- 
tions. We note that this general expectation has already 
been conjectured in Ref. [|l2| on the basis of (n,p) data 
available at that time. These authors have attempted 
to fit the data within a strongly truncated shell model 
calculation which then in turn has been used to pre- 
dict a corresponding (3 decay rate. This procedure is 
viewed as rather uncertain as i) the large energy reso- 
lution in the data made its convolution into a (3 decay 
rate imprecise and ii) the shell model truncation level 
was too inaccurate in order to estimate reliably the con- 
tribution of other states than the back-resonance to the 
decay rate. These shortcomings can be overcome in re- 
cent state-of-the-art large scale shell model calculations. 
We have calculated the (3 decay rates for several nuclei 
under relevant core collapse conditions (p7 = 10 — 1000, 
where pr measures the density in 10 g/cm 3 and tem- 
peratures Tg = 1 — 10). These nuclei include even-even 
ones ( 52 Ti, 54 Cr, 56,58,60^ odd _ A nuclei (59 Mll; 59,6i F6j 

61 > 63 Co) and odd-odd nuclei ( 50 Sc, 54 < 56 Mn, 58 < 60 Co). The 
selection has been made to include those nuclei which 
have been ranked as most important for core collapse 
simulations by Aufderheide et al. ||. In fact, with the 
rates of Ref. || these 15 nuclei contribute between 65% 
and 86% to the change of Y e due to (3 decay in the range 
Y e = 0.44 - 0.47. 

Although the formula for the presupernova [3 decay 
rate is well known (e.g. (|||), we have chosen to quote 
the basic result here as this allows for the easiest discus- 
sion of the improvement incorporated in our calculation 
compared to previous work. Thus, 

f d VV \Q l3 -E e f— ^±ii^) , (i) 
Jo 1 + cxp [kT(p e - E e )\ 
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where E e , p, and [i e are the electron energy, momentum, 
and chemical potential, and C = (Q?f — m^c 4 ) 1 / 2 ; Qif = 
Ei — Ef is the nuclear energy difference between the ini- 
tial and final states, while Sq T is their GT transition 
strength. Z is the charge number of the parent nucleus, G 
is the partition function, G = J2i(2Ji + 1) ex P [ — Ei/kT], 
while F is the Fermi function which accounts for the 
distortion of the electron's wave function due to the 
Coulomb field of the nucleus. The values for the chemical 
potential are taken from fill . 

To estimate the rates at finite temperatures, the com- 
pilations employed the so-called Brink hypothesis PJ24| 
which assumes that the GT strength distribution on ex- 
cited states is the same as for the ground state, only 
shifted by the excitation energy of the state. We have not 
used this approximation, but have performed shell model 
calculations for the individual transitions. Our sum over 
initial states includes i) explicitly the ground state and 
several excited states in the parent nucleus (usually at 
least all levels below 1 MeV excitation energy) and ii) all 
back-resonances which can be reached from the levels in 
the daughter nucleus below 1 MeV excitation energy. As 
these back-resonances also include parent states below 
1 MeV, special care has been taken in avoiding double- 
counting. The partition function is consistenly summed 
over the same initial states. 

Here a word of caution is in order. We have calculated 
the GT strength distributions using 33 Lanczos iterations 
in all allowed angular momentum and isospin channels. 
This is usually sufficient to converge in the states at ex- 
citation energies below E = 3 MeV. At higher excitation 
energies, E > 3 MeV, the calculated GT strengths repre- 
sent centroids of strengths, which in reality are split over 
many states. While this does not introduce uncertainties 
in the summing over the GT strengths (the numerator 
in (1)), it might be inconsistent for the calculation of 
the partition function. However, this is practically not 
the case, as at the rather low temperatures of concern 
here the partition function is given by those states which 
have already converged in our model space. Nevertheless 
there might be states outside of our model space (intruder 
states) which will be missed in our evaluation of the (3 de- 
cay rates. But their statistical weight in both numerator 
and denominator in the rate equation (1) is small. Al- 
though our calculations agree well with the experimental 
informations available (excitation energies and GT tran- 
sition strengths), we have replaced the shell model results 
by data whenever available. 

In Fig. 1 we compare our shell model (3 decay rates 
with those of FFN for selected nuclei representing the 
three distinct classes: even-even ( 54 Cr, 56,60 Fe), odd-A 
( 59 Mn, 57 ' 59 Fe) and odd-odd ( 54 > 56 Mn, 58 Co). We note 
again that (3 decay of the nuclei studied here is impor- 
tant at temperatures Tg < 5 Q. For the odd-odd nuclei 
we calculate rates similar to those of FFN. This approx- 
imate agreement is, however, somewhat fortunate. In 
FFN the misplacement of the GT back-resonances has 



been compensated by too large values for the total GT 
strengths (FFN adopted the unquenched single particle 
estimate) and the low- lying strengths. At higher tem- 
peratures (Tg > 5), the FFN rates for odd-odd nuclei are 
larger than our shell model rates. For odd-A and even- 
even nuclei our shell model rates are significantly smaller 
than the FFN rates as the back resonance, for Tg < 5, is 
less populated thermally than in the FFN parametriza- 
tion. Using the FFN rates, even-even nuclei were found 
to be unimportant for (3 decay in the core collapse; our 
lower rates make them even less important. This sit- 
uation is somewhat different for odd-A nuclei which ( 
57 ' 59 Fe, 59 Mn) have been identified as important in Ref. 
Q adopting the FFN rate. Aufderheide et al. have added 
several odd-A nuclei with masses A > 60 (which are not 
calculated in FFN) to the list of those nuclei which signif- 
icantly change the Y e value during the collapse by (3 de- 
cays. These nuclei include 61 Fe and 61,63 Co; we will show 
below that their rates have also been overestimated sig- 
nificantly in Ref. Q . Our shell model rates indicate that 
the importance of odd-A nuclei is significantly overesti- 
mated when the previously compiled values are adopted. 

In Ref. (l2) f3 decay rates for several nuclei have been 
estimated in strongly truncated shell model calculations, 
in which these authors allowed a maximum of 1 nucleon 
to be excited from the /V/ 2 shell to the rest of the pf- 
shell in the daughter nucleus, and fitted the single parti- 
cle energy spectra to reproduce measured (n,p) data for 
54 > 56 Fe, 58 Ni and 59 Co; the (n,p) data constrain the back- 
resonance transition to the ground states in the (3 decays 
of 54 - 56 Mn, 58 Co and 59 Fe. Our shell model rates are 
compared to the estimates of Ref. jl2| in Fig. 2. For the 
3 odd-odd nuclei, the agreement is usually better than a 
factor 2. This is due to the fact that these rates are dom- 
inated by the back-resonances, i.e. the (n,p) data of the 
daughter nucleus, which are reproduced in our large-scale 
shell model approach and have been fitted in Ref. [fi2"|| . 
For the odd-A nucleus 59 Fe our (3 decay rate is about an 
order of magnitude lower than the estimate of Ref. JlJ] 
at Tg = 2, while the rates agree for Tg > 6, where it is 
dominated by the back-resonances. At the lower temper- 
atures, (3 decays of low-lying states are important which 
might be overestimated in the truncated calculation. 

What might the revised (3 decay rates mean for the 
core collapse? To investigate this question we study the 
change of the electron-to-baryon ratio, Y e , along a stellar 
trajectory. Following Ref. ||, we define 

y » = E Tr** (2) 

k k 

where the sum runs over all nuclear species present in 
the core. Z, A, and X are the charge, mass number, and 
mass fraction of the nucleus, respectively. The mass frac- 
tion is given by nuclear statistical equilibrium B ; we will 
use the values as given in Tables 14-24 of Ref. Noting 
that (3 decay {(3) increases the charge by one unit, while 
electron capture (ec) reduces it by one unit, we have 
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where A^ c and Aj: 



are the electron capture and (3 decay 
rates of nucleus fc. For several key nuclei we have cal- 
culated these rates within large-scale shell model stud- 
ies. Some of the results are listed in Tables 1 and 2, 
where they are also compared to the FFN rates and the 
ones of Ref. §. This comparison also includes the (3 
decay rates for 61 Fe and 61 ' 63 Co, which, due to Q and 
earlier suggested by Brown P,|25|], are important when 
the stellar trajectory reaches electron-to-baryon values 
Y e = 0.44 — 0.46. Our shell model rates agree for 63 Co 
with the rate of Aufderheide et al., but are smaller than 
the estimates of these authors by factors 2 and 5 for 61 Fe 
and 61 Co, respectively. We note that the strong ground 
state decay of 63 Co contributes about 15% to the total 
rate at the condition listed in Table 2. Some of the elec- 
tron capture rates are taken from |14 lq), while several 



other shell model rates are presented here for the first 
time (e.g. for 54 < 56 Fe, 58 Ni and the odd-A nuclei). Al- 
though the nuclei, for which reliable shell model rates 
are now available, include the dominant ones at the var- 
ious stages of the early collapse (due to the ratings in 
Rcf. H), there are upto 250 nuclei present in NSE at 
higher densities ||. Although we are currently working 
at a revised compilation for /3 decay and electron capture 
rates for nuclei in the mass range A = 45 — 65, its com- 
pletion is computer- intensive and tedious. Nevertheless 
some important conclusions can already be drawn from 
the currently available data. 

At first we will follow the stellar trajectory as given in 
Ref. Q, although some comments about this choice are 
given below. We estimated Y^ c and Yf separately on 
the basis of the 25 most important nuclei listed in Ta- 
bles 14-24 in We used shell model rates for the nuclei 
listed in Table 1 and 2. For the other nuclei we scaled the 
FFN rates using the following scheme which corrects for 
the systematic misplacement of the GT centroid and is 
derived by the comparison of FFN and shell model rates 
for the nuclei listed in Tables 1 and 2. The FFN electron 
capture rates have been multiplied by 0.2 (even-even), 0.1 
(odd-A) and 0.04 (odd-odd), while the FFN (3 decay rates 
have been scaled by 0.05 (even-even), 0.025 (odd-A) and 
1.5 (odd-odd). The results for Y^ c ^ are plotted in Fig. 
3, where they are also compared to the values obtained 
for the FFN rates. One observes that the shell model 
rates reduce Y^ c significantly, by more than an order of 
magnitude for Y e < 0.47. This is due to the fact, that, 
except for 56 Ni, all shell model electron capture rates are 
smaller than the recommendations given in the FFN and 
Aufderheide et al. compilations |||5|]. In particular, this 
is drastic for capture on odd-odd nuclei, which due to 
these compilations, dominate Y^ c at densities p-j > 10. 
The shell model (3 decay rates also reduce Yf, however, 
by a smaller amount than for electron capture. This is 
mainly caused by the fact that the shell model (3 decay 



rates of odd-odd nuclei are about the same as the FFN 
rates or even slightly larger, for reasons discussed above. 

It is interesting to note that FFN typically give higher 
/3-decay rates for odd-A nuclei than Aufderheide et al. 
j§, while it is vice versa for odd-odd nuclei. As a conse- 
quence Yf is dominated by odd-A nuclei for Y e < 0.46 if 
the FFN rates are used, while odd-odd nuclei contribute 
significantly if the rates of || are adopted. In cither case, 
both compilations yield rather similar profiles for Y£ C 'P 
(see Fig. 14 in Q). The important feature in Fig. 3 is 
the fact that the j3 decay rates are larger than the elec- 
tron capture rates for Y e = 0.42 — 0.455, which is also 
already true for the FFN rates || . 

So far we have used the same stellar trajectory as in 
Ref. §. This allowed a comparison with the conclusions 
obtained in that reference. However, this assumption is 
inconsistent, and, in fact, was already inconsistent in [^). 
The chosen stellar trajectory is based on runs performed 
with the stellar evolution code KEPLER |^6| which uses 
the FFN electron capture rates, but quite outdated (3 
decay rates ^7j, follo wing the old belief that j3 decay 
rates are unimportant p7[ . The outdated (3 decay rates 
were derived basically from a statistical model approach 
p8| and are orders of magnitude too small j|. What 
are the consequences and will electron capture and (3 de- 
cay rate also balance in a consistent model? At the be- 
ginning of the collapse, electron capture is significantly 
faster than (3 decay (see Fig. 3). The shell model rates 
make 56 Ni the most important contributor, but it can- 
not quite compensate for the reduction of the 55 Co rate. 
Thus, at Y e — 0.485 the total electron capture rate Y" e ec 
drops slightly. This reduction is more severe for smaller 
Y e values, until at Y e — 0.46 electron capture and (3 decay 
balance. The consequence is that, due to the slower elec- 
tron capture, the star radiates less energy away in form 
of neutrinos until Y e — 0.46 is reached. Thus one ex- 
pects that in the early stage the stellar trajectory is, for 
a given density, at a higher temperature. This, of course, 
increases both the (3 decay and electron capture rates. 
Importantly both rates have roughly the same tempera- 
ture dependence in the relevant temperature range: typ- 
ically electron capture rates are enhanced by an order 
of magnitude if the temperature raises from Tg = 4 to 
Tg = 6. But this increase is the same order of magnitude 
by which the (3 decay rates grow in the same temperature 
interval. Consequently the two rates will also be balanced 
at around Y e « 0.46 if a consistent stellar trajectory is 
used. 

As stated above, the dominance of (3 decay over elec- 
tron capture during a certain stage of the core collapse 
of a massive star has been suggested or noted before 
P,p| j2^Jl2] |. However, previous argumentation has been 
based on rates for these two processes which had been 
empirically and intuitively parametrized, rather than de- 
rived from a reliable many-body model. Moreover, it was 
shown in recent years that the assumed parametrizations, 
mainly with respect to the energy of the Gamow- Teller 
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centroid, were systematically incorrect. Shell model cal- 
culations are now at hand which allow, for the first 
time, the reliable calculation of these rates under stel- 
lar conditions. Given the fact that the large-scale shell 
model studies reproduce all important ingredients (spec- 
tra, half-lives, GT strength distributions) very well, the 
shell model rates are rather reliable. We stress an impor- 
tant point, that the shell model (3 decay rates are larger 
than the electron capture rates for Y e « 0.42—0.455. This 
might have important consequences for the core collapse. 
A first investigation into these consequences has been 
performed by Aufderheide et al. |27|, however, using the 
FFN values for both rates. They find that the competi- 
tion of (3 decay and electron capture leads to cooler cores 
and larger Y e values at the formation of the homologuous 
core. These results are important motivation enough to 
derive a complete set of shell model rates and then use 
them in core collapse calculations. 
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TABLE I. Electron capture rates for selected even-even, odd-A and odd-odd nuclei. The chosen stellar conditions reflect 
those at which the nuclei are considered to be most important due to the ranking given by Aufderheide et al. ||. The shell 
model rates (labelled SM) are compared to those recommended by FFN |j| and Ref. The last column, named importance 
ratio, gives the percentage of the total change in Y e (for the definition see text) assigned to the respective nucleus by Aufderheide 
et al. at the respective stellar conditions. All rates are in s _1 . Exponents are given in parentheses. 
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TABLE II. P decay rates for selected even-even, odd-A and odd-odd nuclei. The chosen stellar conditions reflect those at 
which the nuclei are considered to be most important due to the ranking given by Aufderheide et al. [^| . The shell model rates 
(labelled SM) are compared to those recommended by FFN Q and in Ref. The last column, labelled importance ratio, 
gives the percentage of the total change in Y e (for the definition see text) assigned to the respective nucleus by Aufderheide et 
al. at the respective stellar conditions. All rates are in s . Exponents are given in parentheses. FFN did not give rates for 
nuclei with A > 60. 
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FIG. 1. (5 decay rates for several nuclei as a function of temperature and at selected densities at which these nuclei are 
important for electron capture in the presupernova core collapse as suggested by Ref. J5|. The top (middle, bottom) row 
contains even-even (odd-A, odd-odd) nuclei. The solid line shows the present shell model results, the dots give the FFN rates 
II , while the triangles are rates taken from Tables 15-17 in |5J. 
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FIG. 2. Comparison of the present shell model rates for 54,56 Mn, 58 Co and 59 Fe (solid line) with those derived in Ref. 
(dashed line). 
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FIG. 3. Change in the total electron capture and f3 decay rates, Y e ec and Yf, respectively. The shell model results are 
compared with the FFN results along the same stellar trajectory as in Fig. 14 of Ref. pi . 
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